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The critical behavior of the three-dimensional n-vector model in the presence of an external field 
is investigated. Mathematical description is performed with the collective variables (CV) method 
O ' in the framework of the p model approximation at the microscopic level without any adjustable 

^ ! 

^ ' parameters. The recurrence relations of the renormalization group (RG) as functions of the external 

C^ . field and temperature were found. The analytical expression for the free energy of the system for 

-(— > . 
c/2 . 

temperatures T > Tc and different n was obtained. The equation of state of the n-vector model for 

general case of small and large external fields was written. The explicit form of the correspondent 



C^ 



'^ ■ scaling function for different values of the order parameter was derived. The obtained results are 

O . in qualitative agreement with the data of Monte Carlo simulations. 

> 

(N 

en 

o 

in 

o 
o 



X 



I. INTRODUCTION 

The investigation of the critical behavior of the real three-dimensional (3D) magnets is one 
of the most important problems of condensed matter physics. The present work is connected 
with investigation of the classical n- vector model on 3D simple cubic lattice in the presence of 
an external magnetic field by the collective variables (CV) method. Originally this method 
was introduced by Bom^, then used by Zubarev for systems of charged particles^ and later 
developed for calculation of the thermodynamic and structural characteristics of 3D systems 
near the phase transition (PT) point^. The above mentioned model is well-known as the 
classical 0(?2)-vector model or, in field-theoretic language, as the 0(n)-invariant nonlinear 
o"-model. Depending on components of an order parameter this model can describe a number 
of physical systems such as: polymers, ferromagnets, antiferromagnets, the critical point of 
the liquid-vapor transition, the Bose-condensation, phase transitions in binary alloys etc. 

The investigations of the critical properties of the 0(n)- vector models and their partial 
cases were carried out by various methods such as: high- and low-temperature expansions, 
the field theory, the semi-microscopic scaling field theory and Monte-Carlo simulations. 
In general, much attention was devoted to investigation of the universal characteristics 
of the system such as critical exponents and relations of the critical amplitudes of the 
thermodynamic functions. 

The collective variables method as well as Wilson's approach*^ is based on use of the 
hypothesis of scaling invariance and the renormalization group (RG) method for the phase 
transition theory suggested by Patashynskii, Pokrovskii^ and Kadonoflf^. 

The RG method was used to obtain the equation of state of Ising system up to the order 
e^ by Avdeiva and Migdal^ and by Bresin, Wallace, Wilson^. The obtained results were 
generalized for the case of the n- vector model in^. 

Besides, the positive results were achieved at calculation of the thermodynamic functions 
near the critical point. In Wegner's work^- was obtained expression for the free energy with 
taking into account the so-called 'irrelevant' operators in the Wilson's approach. Riedel and 
Wegner— suggested the method of scaling fields for obtaining crossover scaling functions 
of the free energy and the susceptibility. The works of Fisher and Aharony^^, NicoU and 
Albrightfi^ and also Nelson^'* are dedicated to receiving of the crossover scaling functions 
for T > Tc in zero magnetic field near four dimensions. In the frame of the massive field 



theory by Bagnuls and Bervillier— the exphcit results for the correlation length, the sus- 
ceptibility and the heat capacity as functions of the temperature in the disordered phase 
along the critical isochore for one-, two- and three- component systems were obtained. The 
non asymptotic behavior was described as crossover between the Wilson-Fisher's (near the 
critical temperature T^) and mean field's (far from T^) behaviors using three adjustable 
parameters. But this crossover can not realistically describe the situation in the system, 
because there are some physical restrictions of the model. Thus, in the works of Dohm and 
co-workerai^ii^ the calculation of the thermodynamical characteristics of the system without 
e-expansion was performed in the frame of some minimal subtraction scheme based on high- 
ordered perturbation theory and Borel resummation. This minimizing scheme is related to 
the use of the general relations between the heat capacity coefficients for approximation of 
the temperature dependence of coefficient u{t) near the fourth term in the Ginsburg-Landau 
Hamiltonian. This method allows to obtain the non universal critical behavior of the ther- 
modynamic functions below and above the critical temperature, such as the heat capacity 
and the susceptibility as functions of u{t) without any adjustable parameters. But it does 
not provide possibility to analyze the dependence of the thermodynamic variables on the 
microscopic parameters of the interaction potential. 

Besides, the essential success was achieved in calculation of the universal relations of the 
critical amplitudes. Okabe and Ohno^^, Okabe and Ideurap^ investigated the relations of 
the critical amplitudes of the susceptibility by high-temperature-, 1/n- and e-expansion up 
to the order O(e^). Bresin, Le Guillon and Zinn-Justin^^ calculated the universal relations 
of the critical amplitudes for the heat capacity, the susceptibility and the correlation length 
by Wilson-Fisher's e-expansion. 

It should be mentioned, that PT is actively investigated by Monte-Carlo (MC) method. 
Thus Ferrenberg and Landau2i>21 found the critical temperature and critical exponents for 
the Ising and classical Heisenberg models using high-resolution MC method. Besides, the 
universal relations of the critical amplitudes and the equation of state were obtained by 
Engels for 0(1), 0(2), 0(4) models^^""— and by Campostrini et al. for 0(3) model^^. 

A number of new results were obtained using for the description of PT the collective 
variables method. The specificity of the CV method is successive microscopic approach 
and the method of integration of the partition function by short-wave fluctuations without 
applying of the perturbation theory. In the frame of this method the general recurrence 



relations (RR) which correspond to the RG equations, the critical exponents and the relation 
of the critical amplitudes of Ising model were obtained. 

Investigation of the 0{n) model allows to obtain, in the unified form, results for the 
critical behavior of whole class of systems such as: polymers in n — )■ limit, the Ising model 
for n=l, the XY- model for n=2, the Heisenberg model for n=3, the model with n=4 is 
important for quantum chromodynamics with two degenerate light-quark flavours at finite 
temperature and the spherical model in the case n— t- oo which has the exact solution. 

The quantity n is related to dimensionality of the order parameter of the system. The 
investigation of the n-vector model was carried out by the CV method in^I using the 
Stratanovich-Hubbard representation. The CV method was used for investigation of prop- 
erties of the pre-transition behavior and description of the structural PT in the system with 
the n-component order parameter—. The thermodynamical characteristics of the n-vector 
model in zero magnetic field were found in^^''^° using the CV method. 

In general, the real physical systems are characterized by the presence of the external 
fields. The description of systems with the n-component order parameter in the presence of 
the external fields is complicated task and needs detailed study. Thus, taking into account 
the results obtained in2^>^, we investigate the influence of the external field on the critical 
behavior of the n-vector model. 

II. THE MODEL. 

The Hamiltonian H of the n-vector model in the presence of the external field has the 
form: 

i J i 

where Si = {S> , ■ ■ ■ , 5"; ) is the classical n-component spin of length m localized at the N 
sites of d-dimensional cubic lattice with coordinates i, $ (|i — j|) is the interaction potential. 
The partition function of the model ([T]) is the functional integral over all possible orien- 
tations of the spin vector and can be written in the form: m: 

Z = [YI dSiS{Si - m)e-^^, m > 0, (2) 



where we take into account condition that length of the spin is m. We will integrate the 
partition function in the space of CV. Let's introduce the variables : 

^k = ;4E^°^(ki)5i, (3) 

i 



P. = ^^S„ (5) 

i 

which are the n-component vectors. The CV p^^ are introduced as functional representation 
for the operators of the fluctuation of spin density: 

P^ = Jp^Jip-p)idp^f. (6) 

In the CV representation the partition function of the model^i is: 

Z = exp - ^ /3$(A;)pkP-k + ^Po 



J[p]idp^r, (7) 



where 

h = /3H. 

The Jacobian of transition from the spin variables to the CV has the form: 

J[p] = fl[dS-ASi~m)6{po-po)l[S{p'^-pl)S{pi-pi). (8) 

•^ i k 

Calculation of the partition function is performed in the general framework of— i^. The main 
idea is that the phase space is divided on the intervals (layers) to depend on value k and 
the interaction potential is averaged on each of this intervals. The Fourier transform of the 
interaction potential is replaced by the following approximation^"—: 

^^^^ f$(0)(l- 262^:2), keB' 
$(A;) = <^ _ 

[^ $ = const, k G B/B'. 

We assume $ = 0. Such cutting of the potential do not affect the general picture of the 
critical behavior but is appreciable when we want to estimate the critical temperature. We 
use the method suggested in^^ for integration of the partition function in the case of presence 
of the external field. It should be mentioned that in our work we use the quartic measure 



density that allows us to describe the PT on qualitative good leveF^ Integration by 1 layers 
gives: 



where Qi is the partial partition function of I's layer. 



(9) 



go = g^'(w)g^'K), 

Q^\u) = J'[0]expKiV'), 
Qi = Q'''{Pi-i)Q'''{di), 



Q{Pi) = (2 



,(0 



,n + 2 a? 



2 '^^r^piT 



n-i \ fyf 

^,y. exp - 



(10) 



3 (p{xi) 
The function <f{xi) is defined in the appendix. The non integrated part of Z has the form: 



M+i) 



^m 



Yl Pki • • • Pki-^ki. 



.k4 



(11) 



kl...k4<B; + i 

The presence of the external field results in appearance of the linear term in the exponent. 
Let assume that the external field is oriented along one of the coordinate axes (e. g. x axe), 
we receive: 

^m = /(rfPk)^'-exp|v/iV^arV«-i Yl d^'^'\k)p^p_,- 






Y Pki---Pk4'^ki... 



k4 



(12) 



ki...k4<S;+i 

For the coefficients near different powers of Pi^ we have the following recurrence relations: 

C+l) (0 i 



4'+^) = a^^ + d^'\Bi+,,Bi)M{xi), 



,a+i) _ ji).~d 



at>s-~''E{xi), 



(13) 



where 



M{xi) = N{xi) - 1, N{xi) 



yiUrXy. 



XlUnixi 



H. E{xi) 



2dViyi 



^{xi)'' 



Xl 

Ni 



-j^d^^\Bi+i,B), yi = s'Un{xi)J^^^—^, 



(14) 



For convenience the following designation are introduced: 



d^'\Bi+i,B) = rf(')(0) + qs-''\ q = /3$(0)g 



^« --''.,, rf«(0) 



s-'^n, af = s-'W 



(15) 



Thus, the recurrence relations (IT5]1 can be written in the form: 



d+2 

Ui+l = S -i UJi, 



n+i = s'^[{ri + q)N{xi) -q], 
ui+i = s'^~%E{xi). 



(16) 



The initial values of Ui,ri, ui are (for / = 0): 



ojq = s^h , ro = 02 - /3$(0), Uq = a^. 



(17) 



In that way we passed to the parametric space of the RG transformation. The phase tran- 
sition point is represented by the fixed point with coordinates: 



CJ* =0,r* = -/n/5$(0), U* = 0n[/3$(O)]^ 



(18) 



where 



fn = q 



s^[N{x*) - 1] 



s'^N{x*) - 1 
Here x* is the solution of equation^E; 



X 



*2 



N{x*) - s-2 



(19) 



s'^Mv*) = v{x*). 



(20) 



So, when r = 0, h = and / — t- oo the system is in the fixed point. It is obvious that near 
the critical point when r — )• 0, /i — t- and for large 1 in the parametric space the system 



will be near the fixed point. This case is called the critical regime (CR). In the CR the 
recurrence relations may be expanded by the deviation from the fixed point 



n+i - r* 
y M/+1 -u* J 



( *\ 



n 



ri — r 
ui-u* 



(21) 



/ 



The elements of matrix 7^ in linear by {xi — x*) approximation have the form: 

d+2 



Ru 



S 2 



R 



12 



-R22 = VSs^/ii, R 



R2I — Rl3 

„2 



-R3I, 



23 



2Vm* 



{jj,0 - jJ^iX*), 



R32 = V3M*s^-^a;i, i?33 = 
where the following designations were introduced: 



4_d / (^iX 



(22) 



Hi 

tti 
p 



/io Ol 



Qi 



1^0 



n + 2 d 



3(/9(a;*) 



Piy*ri, ri = di 



Qi 



1 



d^^Unjyi) 

dyT 



Univ*: 

ui = ujoih -qi), Uo = s 

1 



y 

2dViy*) 



bi = Qiy*ri, Qr. 



Rfi = i?23V^, R. 



v{y* 



(p{x*) 

<r^v{yi) 

dyT 



(0) 
32 



R. 



32 



(23) 



Action of matrix TZ allows to receive the coefficients of the partition function of the next 
layer. In order to receive the coefficient of I's layer we need to act by TZ 1 times on the 
coefficients of zero layer. 






ri — r 
ui-u* 



^ LUo-OJ*^ 



7^' 



/ 



To — r 
Uq — u* 



(24) 



/ 



It is easy to obtain the form of the matrix 7?.' if reduces the matrix 7^ to the diagonal form. 
In order to do it, we must pass to the base from the eigenvectors of 7^. The eigenvalues of 
TZ are universal quantities: 

1 



El = Rii, E: 



11; -C/2,3 



-R22 + -R33 ± V (-R22 — -^33)^ + 4i?23-R: 



32 



(25) 



The eigenvectors have the form: 



Wi 





VV 



^2 



/o \ 

1 



U-s 



R 

VJ 



(26) 



The inverse vectors are written as: 

Vi = (10 0), V2 



D 



(0 1 -R), vs = -{0 -Ril) 



D 



R 



R 



23 



Ri 



Eo — R 



22 



The determinant of the inverse matrix is 

D = 



R 



23 



E-i — Eo 



E3 — R22 
After expanding the coefficients in (IT6l) by eigenvectors we obtain: 

- ; 

cji = s^hEl, 

ri = r* + ciEl + C2REI, 

ui = u* + CiRiEl + C2EI 

The coefficients ci, C2 can be found from the initial conditions at 1=0. 

1 



(27) 



(28) 



(29) 



Cl 



D 

1 



[ro-r* - R{uq -u*)], 



C2 = -^[{uq - u*) - Ri{rQ - r*)]. 



(30) 



The eigenvalue i?2 > 1 is responsible for the deviation from the fixed point, and i?3 < 1 is 
much more smaller than E2 and we can neglect it. This approximation neglects the confluent 
corrections. Taking into account that at the phase transition point ri = r*, we obtain the 
equation for the critical temperature in the forn>^: 

[/3,$(0)]2(1 - /„ - R*y^) - a2/3,$(0) + a4i?7v^ = 0, (31) 

R* = Ry^. (32) 

This equation allows to write the solutions of the recurrence relations as function of the 
temperature and the external fleld for the critical regime: 



ri 

Ul 



4h'E[, 



/3<l>(0) (-/„ + CittE^ + R*C2TEy^/^ 
[/3<l>(0)]' (0„ + c,Tr^/^RlE', + C2tE[) . 



(33) 



The obtained coefficients are exponential functions of 1. For small 1 their values are small 
and then increase rapidly. For big 1, the value of ri (the coefficient near the square term in 
the partition function) is bigger than Ui (the coefficient near the quartic term in the partition 
function). Thus, we can integrate the partition function in this region using the Gaussian 
approximation. But in the CR we must use the distributions of fluctuations higher than 
Gaussian's one and take into account the quartic term in the partition function. So let's find 
the number 1 after what we can pass from accounting quartic to accounting only quadratic 
terms. We call this number as the exit point from the CR. In zero magnetic field this point 
was studied in^. Let designate it as m,- and find it from the condition of deviation from the 
fixed point: 

rm.+i-r* = ~6r*, 6 = ^ (34) 



As a result we obtain: 



where 



'* = -5r*, 6 = 


T 




T 


In f 
hiE2 




Clk 

T = T 

Jn 





"^- = -^T-ET - 1. (35) 



(36) 



is the renormalized reduced temperature. In the presence of the external field the exit point 
from the CR is found from the condition: 



From this we found^: 



Wn^+l -u* = ho. (37) 



\nh - !ih' . , 

n. = -j^-l, /. = .o^. (38) 



The quantity h^ is found from the normalization condition for scaling function. 

If there are both the temperature and the field the exit point depends on relation between 
T and h. Some boundary temperature field he which divides the values of fields on strong 
and weak was found in^*^. The condition of this division on strong and week fields has the 
form: 

rrir = nh. (39) 

After substituting equations for the exit points from the CR we obtain: 

K = |r|^°, (40) 



10 



where the critical exponent pq has the form: 



InEi u 2 



The critical exponent of the correlation length u ioT h = is 

Ins 



z/ 



The critical exponent /i describes the dependence of the correlation length from the field for 
T = Tc. Therefore we can rewrite fl35|) in the form 

m. = -^-^ - 1. (43) 

In general case deviations by the temperature and by the field from the fixed point should 
be united. This united point Up can be found from the equation^: 

(sJ/'/iX'+')' + (ciTrmO)E2'^'y = r*^- (44) 

The value of Up can be found numerically. But numerical solution does not allow analyti- 
cally to take into account influence of the temperature and the field on the critical behavior. 
Thus, in^^ formula for united exit point from the CR was introduced. This form includes 
the temperature and the field and allows in the limit of zero field or zero temperature pass 

to the Eq.([38]) or Eq.g3]) 

In {P + hi' 



Ci 



Information about the exit point from the CR allows us to divide the integration of the 
partition function on two stages: integration by quartic distribution in the CR, where values 
of Ti and ui are commensurable and integration by Gaussian distribution (the Gaussian 
regime) for the phase space layers with I > np. But deviation from the fixed point occurs 
not so sharp to exactly distinguish the critical and the Gaussian regimes. Therefore we 
should take into account the transition regime (TR) in which r^ exceeds Ui but we still can 
not use Gaussian distribution. In order to integrate the partition function in the TR one 
must use the quartic measure density. Fortunately, there is only one layer in the phase space 
for which r; and ui behaves like described for the TR. The number of this layer is next after 
the CR Up + 1, so from the layer with Up + 2 the Gaussian region is present. Thus the 
partition function has the form: 

Z = QoQi ...Qn, Qn,+i 2^('^"-+^-')g'^"''+HP„,+i)^„,+2, (46) 

CR TR LGR 

11 



where CR designates the critical region, TR is the transition region and LGR is the hniiting 
Gaussian region. 

The coefficients near the CV change their behavior for / > n^. It siniphfies calculations. 
Let's take / = rip + 1 in flTU|) . Then in order to integrate Z„ _|_2 is important to know below 
or above T^ is the system. The coefficients 






(47) 



depend on the temperate and the field. The values of m„ _|_2 are always positive that provide 
convergence of (H6|) . The coefficient r„ +2 is positive and exceeds m„ _|_2 for large r {h^ » 
h). In this case the partition function can be integrated in Gaussian approximation. But 
for small r [h^ « h), rnp+2 decreases and becomes negative. In this case the Gaussian 
approximation is useless. This problem can be solved by introduction the substitution: 



Pk = < + ^<^k, a+ = (a+,...,0). 



(48) 



Then ^n„+2 changes to: 



where 



Zn,+2 = e^^°('^+) f{dr]f"p+^ exp 



i J2 d^''^^'\k)rj^rt_^- 



keBi 



np+2 



(np+2) 
ri 3 3(nj,+2) 

12 ° 

(np+2) 
(np+2) 



^' Yl ^k^-k + 2 J] r^i'^r/ 



(i)Ji) 

k 



5 2gf(np+2)^^ 



5Z ^& 

kl...k3GBrip+2 



^k2%3^ki+k2+k3 



24K 



np+2 



Yl ^ki • • • ^k4'^ki+...+k4 



ki...k4eB„p+2 



04 



Eo(^+) = h'a^ - ;-d"^^\0)al - ^sls'^-^^'^a 



24 



The shift cr+ is found from the condition: 

dEoia+) 



da 



0. 



+ 



(49) 



(50) 



(51) 
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This condition causes the coefficient near the first power of tJq equals to zero and results 
to a cubic equation for cr+. Solution of this equation can be found in the form 

a+ = aos("''+2)/2^ (52) 

Finally one obtains: 

cTl+p(7o + q = 0, (53) 

where the following designation were introduced: 

j,.^^,,.__^_^. (54) 



Solutions of a cubic equation depend on the sign of the discriminant 



^ (IJ H- (Ij . (55) 



For T > Tc the value of Q is always positive. So the Eq. (153|) has one real and two complex 
roots. We select the real one: 

ao = A + B, 

A={-l+VQy.B={-l-^)\ (56) 

The integral (H9|) is calculated in Gaussian approximation. This integration terminates the 
calculation of the partition function: 






loo 



(n-l) 

di{k) \d2{k) 



1 / 1 \ \"~^^ 



where: 



d,{k) = rf + (5^{{])h''k\ 
rR - s 2 ' 

~(2) 
^(2) _ „-2(n„+2)^fl 

'^ii - 5 2 ' 

4^^ = '^np+S + 2^0^0Wnp+2, 



r 



^R^ = ^np+2 + gS|^a^M„p+2. (58) 
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Thus, dividing the phase space of the CV on layers depending on the values of the wave 
vector, we distinguished two main states of the system to refer to the critical behavior. 
The critical region which corresponds to short-wave fluctuations and the Gaussian region 
which corresponds to long-wave fluctuations. Short-wave fluctuations are connected with 
the microscopic parameters of the model and long-wave fluctuations determines the critical 
behavior. In contrast to earlier approaches we pay equal attention on both 'irrelevant short- 
wave' and long-wave fluctuations. 

III. THE FREE ENERGY 

After calculation of the partition function we can flnd the free energy of the system. 

F = -kT\nZ. (59) 

Described above structure of the partition function allows to present the free energy as a 
sum of terms which corresponds to different regimes of fluctuations: 

F = Fq + FcR + Ftr + Flgr- (60) 

Here 



,^,,, ■(27r)tm"-i 
Fo = -kTN\n ^ — ; , , 

T{n/2) 



(61) 



is the free energy of noninteracting spins, 

rip 

FcR = -kTj2^nQ, (62) 

1=0 

is the energy to refer to the critical region, 

Ftr = -kTQ^^+^ (63) 

is the contribution which corresponds to the transition region and, respectively 

Flgr = -kT\n [2t(^"p+i-^".)Q^"p+2(P„^+i)^n,+2] (64) 

is the contribution which comes from the Gaussian region. 

Eq. fl62p depends on the exit point from the CR. In order to know this dependence explic- 
itly we need to sum up elements Qi. For this purpose let's extract dependence on index 1 in 

14 



Qi. The quantity Qi is a function of t//. The variable y is bigger 1 {yi » 1) for any tem- 
perature. Thus, we can use the expansion for Weber's parabohc cyhnder function U{a, x) 
in Qi by inverse powers of yi. Using appropriate expansions we have: 



FcR = -kTN'f^^-kTj2Nifi: 



(65) 



/=i 



where 






(2n + 7), 



.'-' + 



71 — 1 \ 71 



n — I 



ln( ^l+l^f- 



(66) 



To extract the dependence on / from fi we expand it by powers of Xi — x* and substitute in 
obtained formula the solutions of the recurrence relations. After summing up we obtain: 



CR 



-kTN' 



ioi+li^ + l2r'-l\h' + KYo 



where 



i - 


= 71 + 72 iJc + 73^c , 


71 = 


f* 

J cr 
l-S-3' 


72 = 


fndi5^ 


73 = 


PndsS^ 



For the free energy of the TR we obtain: 



(67) 



(68) 



For the GR we have: 
InZ 



TR 



-kTN'U^^,{h' + hl)l., 



Jrip 



+1 



Inf/ 



n — 1 



• ^rip+l I + 



X. 



■rip + l 



+ 



np+2 






71 TT 1 TT 



n — 1 TT 

H In 



Ei-Z7I:T + ("-1)Ei-Z7 



4(0) ^ d,{k) 



fc>0 



d2{k)' 



(69) 



(70) 
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In order to sum up by k we change the summation by integration, that results: 



where 



'LGR 



^(+) 



kTNr 



p+2 <^ n 



1 



In 2 + In s In 3 + - In -u 

4 4 



rip+l' 



1, ... . n + 2 



1 

2 L 



lnr2^ + /^, + (n-l)(lnrt) + /^^ 



/A,, = ln(a^ + 1) 1 — K T arctan Oj, 

3 fli fli 



7r6 //3$(0) 



Co 



f(0 • 

' R 



For the convenience the following designations were introduced: 



(71) 



(72) 



LGR 



4^' + Fa, 



(73) 



where 



F^-^^ = -kTNEoia^) 



Fg = -kTN^^^^fc, 
fc = n 



1 1 1 

-ln2 + lns- -ln3 + -\nUnp+r 



1, ... ^ n + 2 
-lnf/(x„^+i)-— 2 

^ ^Vnp+l 



1 



^(1) 



lnrt^ + /G. + (n-l) lnrt^ + /^, 



^(2) 



(74) 



Finally, we obtain the free energy which is dependent on the temperature and the field: 



where 



-kTN I In 



(27r) 2 171 



n— 1 



:(7oi + 7i^ + 72r^)- 



- eoh' (h' + hi) "' - (7+ - 62) (h' + h 



(75) 



eo 



62 



f^o 



0"n 



2s3 



p+2 + T^'"np+2So' 



v^ 



7; 



i (^"^^1 - ^' + §) • 



(76) 
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IV. THE ORDER PARAMETER. 

From the formula for the free energy we obtain the order parameter of the system by 
direct differentiation by the field: 

The structure of the free energy allows to separately differentiate parts connected with 
different fluctuation processes. The result reduces to the form: 

M = aUh' + hl)^o. (77) 

The quantity ct^Jq depends on the variable a which represents the ratio between the field and 
the temperature: 

1 a'^ \ a 



a = h/K. (78) 

Below there are equations for the coefficients that depend on a, n and other parameters 
of the model. The method of calculation of the order parameter for the n-vector model is 
analogous to the method for the Ising niodel'^^. For convenience we use the same designations, 
but in n-vector model appears dependence on n: 

eo = -^, 

R 3/2 



eo2 = ^ /71 + f^o?. 



1 + Y^^A 
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cd 



4 1 + $ffcd 



+ 



4 <+i 2U{Xn,+l) 



9p+i 



C + <^^« + (--l)fS + af^f 



.(1) 



t(2) 



7p 






n 
2 






HU72 + 2%H,). 



(79) 



Thus we obtained explicitly the equation of state of the n-vector model in the presence 
of the external field. Its form allows easily to pass to boundary cases of dependence only 
on the temperature or on the field. So this equation is called the crossover equation. The 
quantity ct^Jq is the scaling function of the crossover equation of state. It depends on the 
ratio of the field to the temperature a. Eq.( 1771) allows to obtain a graph of dependence of 
the order parameter on the field for T = Tc and compare it with results of Monte-Carlo 
simulations for analogous models. Graphs of such dependencies for different n and different 
parameters of the interaction potential are presented on figiH El |3l HI As we can see from 
these figures the order parameter M decrees when the componence of the model n increases. 

There are different forms of the equation of state. Some discussion about convenience the 
correspondent forms of the equation of state is presented in^. The equation of state ( JTTll 
can be reduced to the form used in^^i""— : 



M = h'/'fciz), 



where: 



h = H/Ho, z 



t 



t 



T. 



(80) 



(81) 



/i/35' " ' To 

and fc is the scaling function. The explicit form of fc can be found from extrapolation 
of the MC data obtained idM'^—. Hq and Tq are normalization constants. Such form is 
equivalent to the Widom- Griffiths equation of stated 



y = f{x) 



(82) 
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Figure 1: The dependence of the order parameter on the field for r = and n=l 
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Figure 2: The dependence of the order parameter on the field for r = and n=2 



where 



y = h/M\ X = t/M^I^. 
The scaling function obtained from Eg. (1771) has the form 



(83) 
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Figure 3: The dependence of the order parameter on the field for r = and n=3 
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Figure 4: The dependence of the order parameter on the field for r = and n=4 



It depends on a 



a 



h 

fPO ' 



(85) 



The variables a and z are connected with the ratio: 



a 



„3/2 / r \po 

^0 I Jn \ ^-po 



ho V Cifc 



(86) 
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Figure 5: The scaling function for n=l and h/c = 0.5, solid curve - our results, dashed curve - 
Monte-Carlo data;^^. 

that allows to compare our results with Monte-Carlo data. Figures [5l|8] present graphs of 
the scaling functions for different n, where the dashed curve is the Monte-Carlo data. 

V. CONCLUSIONS. 

We obtained the partition function of the n- vector model in the presence of the external 
field above the critical temperature by the CV method. The method of calculation corre- 
sponds to the general scheme of the RG approach. Taking into account the explicit form of 
the interaction potential allows to obtain the explicit dependence of the coefficients of the 
linearized recurrence relations on the temperature and the microscopic parameters of the 
model. 

The explicit form of the exit point from the CR allows to obtain the equations for the 
recurrence relations in the CR suitable for any ratios of the temperature and the field. The 
structure of the partition function as a product of the partial partition functions that present 
different fluctuation processes allows to obtain the explicit form for the free energy of the 
system. The order parameter of the model was found by direct differentiation by the fleld. 

The formulas to describe the fleld dependencies of the order parameter of the n-vector 
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Figure 6: The scaling function for n=2 and h/c = 0.5, solid curve - our results, dashed curve - 
Monte-Carlo datai^i. 

model with exponentially decreasing interaction potential for different ratios h/c {h is range 
of the interaction potential, c is period of the simple cubic lattice) were obtained. It was 
found that for each value of n there is an appropriate value of h for which our dependencies 
are close to Monte-Carlo data (see FigJSHH]). 

The explicit form of the scaling function (18^ was found. The comparison with Monte- 
Carlo data shows some difference for the behavior of fciz) for intermediate values of z. It 
may be caused by used approximation in which the critical exponent 77 = and corrections 
to scaling were neglected. But further specification of calculations is a subject of a separate 
investigation. 



VI. APPENDIX 
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Figure 7: The scaling function for n=3 and h/c = 0.5, solid curve - our results, dashed curve 
Monte-Carlo data^. 
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Figure 8: The scaling function for n=4 and h/c = 0.5, solid curve - our results, dashed curve 
Monte-Carlo data2^. 
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U{a, x) is the function of Weber's parabolic cylinder. 

— 

3 \ ^^ fn — 1 \ /x 



Qidi) = (27r)- \—^\U[ ^—,xi exp 



(0 / ^ V 2 '-^'/^^^^ V 4 



a4 = -3s^— (1 - z'Uniz') - U^) , 



m"^ 



Uo = Xr^UrXz'), Z' = ./4«2, 

V 2 y M4 



' p+m 



9p+m 




xE^-^H,d L H.^El^-^ 



1 + <^E^-'H, \ 2 [1 + $E2"^-iiJj 



X = fn 
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